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SUMMARY 


This report provides an updated version of analytical documentation for the V072 
Rotor Wake/Stator Interaction Code. It presents the theoretical derivation of the equations 
used in the code and, where necessary, it documents the enhancements and changes made 
to the original code since its first release. 

V072 is a package of FORTRAN computer programs which calculate the in-duct acoustic 
modes excited by a fan/stator stage operating in a subsonic mean flow. Sound is generated 
by the stator vanes interacting with the mean wakes of the rotor blades. In this updated 
version, only the tonal noise produced at the blade passing frequency and its harmonics, is 
described. The broadband noise component analysis, which was part of the original report, 
is not included here. The code provides outputs of modal pressure and power amplitudes 
generated by the rotor- wake/stator interaction. 

The rotor/stator stage is modeled as an ensemble of blades and vanes of zero camber 
and thickness enclosed within an infinite hard-walled annular duct. The acoustic pressure 
within the duct is calculated by distributing pressure dipoles on the surface of the stator 
vanes and calculating the pressure at an arbitrary point within the duct via the normal 
mode expansion of the Green’s function for an annular duct. By this procedure one obtains 
an infinite series for the sound pressure within the duct. Each term contains a normal mode 
of the duct multiplied by the amplitude of that mode. The amplitude of each propagating 
mode is computed and summed, with appropriate factors, to obtain the harmonics of sound 
power flux within the duct. These calculations are carried through for both upstream and 
downstream propagating modes. 

No assumptions regarding the ratio of the wavelength of the sound generated to the 
vane chord are required. However, to simplify the computation of the dipole distribution 
generated by a given periodic inflow variation, the so-called strip theory approximation is 
employed. Specifically, the dipole distribution is computed by deleting, from the convected 
wave equation, all terms containing derivatives with respect to radius, but retaining the 
radius as a parameter in the boundary conditions (which derive from the vane geometry 
and the incident fluid flow). At each radius, therefore, the equations to be solved are 
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those for a linear cascade in subsonic flow. Although radial derivatives are deleted from 
the wave equation in calculating the chordwise pressure distribution on the vanes, radial 
variations in the amplitude and phase of the pressure distribution are taken into account in 
integrating over the vane surface to obtain the amplitudes of the propagating duct modes. 

The wake flow downstream of the rotor is modeled as a small amplitude disturbance 
flow, superimposed upon a steady mean flow. The mean part of the flow is assumed to 
have no radial component, so that the mean flow stream surfaces are cylindrical. It is 
further assumed that if one of these cylinders is unwrapped to form a plane, the mean flow 
streamlines will be parallel. This is equivalent to ignoring viscous diffusion in the rotor 
wakes, a justifiable approximation as regards the calculation of the forces on the stator 
vanes, provided that the wake thickness chosen is that which is obtained at the axial station 
of the stator vanes themselves. The magnitude of the wake flow velocity is constant on 
linps parallel to the mean flow streamlines, but varies periodically in the direction normal 
to the mean flow streamlines. 
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CHAPTER 1 


INTRODUCTION 


1.1 Background 

This report is a revision of the original analytical documentation for a rotor-wake/ sta- 
tor interaction computer code described in NASA contractor report, CR-167952. The code, 
which was originally developed under contract to NASA by Bolt Beranek and Newman 
Inc., was later revised and integrated into a fan noise prediction system by Pratt &; Whit- 
ney and renamed the V072 code. This report replaces Volume 1 of the original two- volume 
documentation authored by Ventres et al. 1 Volume 2, “Computer Programs,” of the same 
set 2 is superseded by a revised technical documentation and user’s manual written by 
Topol and Mathews. 3 The current revision, while primarily a mathematical treatment, 
should also enable a diligent reader to trace through the actual computer code without the 
need for the old Volume 2. 

This revision was necessary for several reasons. First, and most important, because 
the original documentation preceded the actual release of the computer program, it did 
not reflect many of the changes to the code. In particular, the geometry used for the code, 
and, as a consequence, the resulting coded formulae, are different from those described in 
the old Volume 1. Second, the original report covered, in addition to discrete-tone noise 
sources produced by mean rotor wakes, broadband noise sources due to turbulent rotor 
wakes. These broadband noise capabilities are included in other codes but are not part of 
the V072 computer program. Finally, there have been several undocumented corrections 
and enhancements to the code since it was first released. 

Outside of the changes described above, most of the material in this revision has been 
borrowed intact from the original documentation, although some of it has been rearranged 
to improve the flow of the text. Where necessary, complete sections have been deleted 
from the material. These include all of the discussion on broadband noise sources which, 
as mentioned earlier, is not included in the computer program. These sections cover 
mode lin g of the rotor turbulent wakes and their attendant acoustic response computations 
(i.e., broadband noise spectrum and acoustic power). Therefore, with no further need for 
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a general frequency treatment, all of the material in this revision is rewritten in terms of 
discrete harmonics of the rotor blade passing frequency. 

1.2 Overview 

Given a description of the rotor and stator geometries and their operating conditions, 
the V072 code calculates the modal content of the acoustic fields set up inside the inlet 
and exhaust ducts of a turbofan by the actions of unsteady vortical velocity fluctuations 
convected past the stator vanes. The results from this code, supplemented by mathematical 
models of the impedance of the inlet and exhaust terminations, or coupled with stand- 
alone radiation codes such as the one developed by Eversman , 4 can be used to calculate 
the farfield noise. The program can also be used to provide the initial conditions required 
for the investigation of duct liners, particularly liners tailored to suppress a specific mode 
or set of modes. 

In the context of discrete rotor/stator interaction noise, the fluctuating velocity field, 
i.e., vortical flow, is caused by the periodic wakes of the upstream blade row of a fan 
operating in front of an outlet guide vane. The fluctuations are assumed to be small 
enough so as to be convected by the mean flow. When these velocity fluctuations are 
carried past the rigid surface of a stator vane, the requirement of flow tangency sets up 
acoustic pressure fluctuations, which propagate upstream and downstream in the duct, 
and are perceived as noise. 

The acoustic pressure in the duct can be calculated in a straightforward manner using 
the Green’s function for the duct, provided that the pressure distributions on the surface 
of the stator vanes are known. If the Green’s function is expressed as an infinite series 
of the normal modes of the duct, the modal amplitudes of the acoustic pressure field 
are obtained directly. The problem remaining, therefore, is to determine the pressure 
distributions on the surface of the stator vanes generated by the convected vortical flow. 
This is accomplished by modeling the vanes as surfaces of zero thickness and camber, 
which support an u nkn own but continuous distribution of pressure dipoles. The normal 
velocity on the vanes induced by the unknown dipole distribution is then required to 
nullify the normal component of the vortical disturbance velocity so that the total fluid 
velocity conforms to the shape of the blades. The result is an integral equation for the 
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dipole distribution, which, depending upon the approximations introduced, may admit an 
analytic solution, or may require numerical treatment. 

The appro xim ations required to arrive at an analytic solution are considerable; es- 
sentially they ignore the annular or circular geometry of the fan, or stator, the three- 
dimensionality of the disturbance, the interaction between the blades, and the compress- 
ibility of the fluid. What remains is incompressible flow around an isolated airfoil in two 
dimensions, the solution of which, for a sinusoidal vortical gust, is known as the Sears func- 
tion. This was the approach used by Kemp and Sears 5 ’ 6 in their original investigations 
of the interactions between blade rows in axial flow turbomachinery. The incompressible 
approximation has since been shown to be unacceptable. For example, Fleeter 7 found that 
compressibility can change the pressure values on the rotor or stator by as much as a factor 
of two, while Kaji 8 calculated even greater changes (20 dB; a factor of 10) in the sound 
pressure level upstream of the blade row due to source non-compactness. 

At the opposite extreme, as regards complexity, are calculations carried out by Koba- 
yashi 9 and Kobayashi and Groeneweg 10 , using equations derived by Namba 11,12 for an 
annular blade row in compressible flow. Because of the annular geometry, the kernel 
function, which relates the dipole distribution to the velocity normal to the blades, is elab- 
orate in form and time-consuming to compute. The associated integral equation is two- 
dimensional (i.e., both the spanwise and chordwise distributions of pressure dipoles must 
be determined). In his numerical work, Kobayashi made comparisons between “exact” 
aerodynamic theory and various approximate methods of calculating the normal compo- 
nent of the induced velocity. One of these, which he called the “quasi-three-dimensional” 
approximation, coincides with what most would call a “strip theory approximation. In 
calculating the induced velocity, derivatives with respect to radius are deleted from the 
convected wave equation, so that one is left with the kernel function of a linear cascade in 
unsteady compressible flow. The particular geometry of the cascade, as well as the vor- 
tical inflow velocity, depend on the radius, and this parametric dependence is retained in 
calculating the induced velocity and, hence, ultimately the dipole distribution. The effects 
of radial or spanwise variations in the amplitude and phase of the vortical disturbance 
velocity, as well as radial variations in blade chord, inter-blade gap, and stagger angle, 
are thereby taken into account, even if in an approximate manner. Kobayashi found that 
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this procedure introduced errors no greater than 2 dB in the computed magnitudes of the 
acoustic modes set up in the duct upstream and downstream of the rotor. Furthermore 
the strip theory approximation, as described above, is considerably easier than the “exact” 
method to implement numerically. For these reasons, this approximation was selected for 
use in the V072 program. 

Laksminarayana et al. 13-15 have collected an extensive set of measurements of the 
mean and fluctuating components of the flow downstream of a multi-bladed fan operating 
at subsonic tip speed, using a transducer rotating at the same rate as the fan. Their 
measurements of the mean, or time-averaged, component of the flow at locations remote 
from the blade hub or tip indicate that the radial velocity immediately downstream of 
the fan blade trailing edges is substantial, but that this velocity decays rapidly, so that 
at distances greater than about one-half blade chord downstream the circumferential and 
axial components of the mean velocity predominate. Measurements of the velocity defect 
profile for each fan blade show marked asymmetry about the streamlines on which the 
minimum velocity occurs, but this asymmetry also disappears within a half-chord length 
downstream of the fan blade trailing edges. At greater distances, the normalized velocity 
defect profiles all show a Gaussian distribution.* Due to the diffusion of momentum in the 
flow, the widths of the velocity defect profiles increase with downstream distance, while the 
ma ximum velocity defect decreases. But for the purpose of calculating the noise radiated 
by the stator, it is sufficient to model the wake flow in the vicinity of the stator. This can 
be done by ignoring momentum diffusion between the fan and stator, provided that the 
wake thickness and velocity defect are assigned the values they attain in the vicinity of the 
stator, say at the leading edges of the stator vanes. 

With these approximations, the mean flow streamlines downstream of the fan become 
a series of parallel lines, as viewed on a cylindrical surface opened out to form a plane. 
On this plane the velocity defect profile of each rotor blade is assumed to be a Gaussian 
wake profile, with wake width and maximum velocity defect left as parameters to be 

* Since the original version of the code was released, hyperbolic 16 and loaded rotor 17 
wake profiles, which are more representative of modern fans, have also been added. The 
loaded rotor model, in particular, indicates that asymmetric profiles may have a non- 
negligible influence on the predicted noise levels. 
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specified. When viewed in a stationary frame of reference fixed to the stator, the pattern, 
which is spatially periodic, becomes a temporally periodic fluctuation which can most 
conveniently be expressed in terms of a Fourier series for the purpose of calculating the 
inflow disturbances normal to the stator vanes. Further details of the implementation of 
the wake model in the code will be given in Chapter 3. 

One additional aspect of the mean wake geometry is noteworthy. At each radius the 
mean flow streamlines all bear the same angle to the centerline of the duct, but this angle 
is a function of the radius, the precise nature of which depends upon the radial variation 
of the rotor loading. The locus of centerlines of any given wake form a surface whose shape 
depends upon the radial variation of the wake angle; this surface will normally intersect 
any selected stator vane at only one point, and this point will move along the leading edge 
of the stator vane as the fan rotates. If the gap between the fan and stator is large, each 
stator vane will at any given instant intersect many rotor blade wakes, and, for each wake 
intersected, there will be two changes in the sign of the disturbance velocity seen by the 
stator vane. The pressure induced on the stator vane will have as many sign reversals as 
the inflow velocity, and will tend to excite duct acoustic modes with that same number of 
radial nodes. If, at the frequency of interest (which, of course, must be the blade passage 
frequency or one of its harmonics), few or none of the duct modes with that number of 
radial nodes propagate in the duct, the stator will not, at that frequency, be an effective 
generator of sound. The significance of this wake roll-up phenomenon in the generation 
of sound by the stator was pointed out previously by Bliss, et al. 18 In their discussion, 
they chose to emphasize the radial trace velocity of the vane/ wake intercept instead of 
the number of nodes in the normal component of the inflow velocity. Their criterion for 
efficient sound generation, namely that the trace velocity of the points of intersection of 
the wake centerlines with the stator vanes be supersonic, is perhaps more appropriate to 
an unshrouded rotor/stator combination, but the two criteria are roughly equivalent when 
the ratio of hub and tip radii is close to one, and are exactly equivalent for a set of two 
lin ear cascades in relative motion between two infinite parallel planes. 

In the equations derived in this report, and in the computer program developed from 
them, the radial variation of the angle between the mean wake streamlines and the duct 
centerline may be specified arbitrarily, so that the effects of wake roll-up are properly taken 
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into account. 


One element of the program, which was absent in the original version, is a model for 
rotor hub and tip vortex flowfields. The details of the model may be found in Ref. 16, 
where measured rotor flow data were used to develop empirical relations for the vortex 
flowfield. While the model may not be sufficiently robust to allow for realistic calculations 
of the vortex noise contribution, it does enable one to perform a parametric assessment of 
the potential impact of a tip, or hub, vortex on the rotor/stator interaction noise. 

The remainder of the report is organized as follows. Pertinent aspects of the duct, 
rotor and stator geometries are discussed in Chapter 2. Mathematical models of the fan 
blade mean wake and its hub and tip vortex flows, which together represent the rotor 
gust upwash, are presented in Chapter 3. The chapter then concludes with a discussion 
of the computation of the stator vane pressure distribution generated in response to the 
gust upwash. An overview of the acoustics of annular ducts is given in Chapter 4, along 
with a derivation of the amplitudes of the duct modes excited by stator vane pressure 
distributions. The expression for the acoustic power in the duct is also presented. Chapter 
5 presents key formulae for duct mode amplitudes and acoustic power derived in earlier 
chapters in non-dimensional forms corresponding to those used in the code. Details on 
cascade response and numerical methods are given in Appendices A through E. 
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CHAPTER 2 


GEOMETRY OF DUCT, ROTOR AND STATOR 

For our purposes, a turbofan is a rotor and stator combination mounted in an annular 
duct of infinite length (i.e., reflection of acoustic waves from the ends of the duct are 
ignored). A sketch of the rotor/stator combination is shown in Fig. 2.1. As shown in the 
figure, a cylindrical polar coordinate system is established in the duct, with the polar axis 
lying along the duct centerline. The axial coordinate Xi, increases in the direction of the 
air flow, and the rotor rotates in the sense of increasing 4> with fan rotational speed Cl. 
Note that the meaning for X]., X2 and X3 in subsequent figures is sometimes different from 
that in Fig. 2.1. The outer radius of the duct is r D and the inner radius is r H . The rotor 
has B identical evenly spaced blades and the stator has V identical evenly spaced vanes. 
This is the same coordinate setup as used in Ref. 19, page 190, from which we later obtain 
our Green’s function. 

The rotor blades and stator vanes are modeled as twisted sheets of zero thickness and 
camber, whose stagger and chord vary with the radius, r. The axial and azimuthal sweep, 
if any, of the stator vanes are defined by extending a radial line from the axis of rotation 
through the leading edge point of the vane at the hub (radius r = r H ). As shown in Fig. 
2.2, the location of the leading edge point at any other radius r is defined by the two 
parameters x SD and y SD which give the displacements of this point from the radial line. 
(x SD = axial displacement of the stator vane leading edge, y SD = azimuthal displacement 
of the stator vane leading edge.) Similarly, the azimuthal sweep of the rotor blades is 
defined by extending a radial line from the axis of rotation through the trailing edge point 
of the blade at the hub. The location of the trailing edge point at any other radius r is 
defined by the parameter y RD which gives the displacement of this point from the radial 
line. ( y RD = azimuthal displacement of the rotor trailing edge.*) In general, x SD , y SD and 
y RD are functions of the radius, and by definition, x SD = y SD = y RD = 0 at the hub. 

* No axial sweep parameter is used for the rotor. Instead a rotor/stator distance pa- 
rameter, x SPAC . which represents the axial distance between the rotor trailing edge and 
the stator leading edge, accounts for any axial sweep that might exist. 
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Fig. 2.2 Blade and Vane Geometry. 




The intersections of the rotor blades and stator vanes with a cylindrical surface of 
radius r is shown in Fig. 2.3. This is for the point in time where the rotor trailing edge 
and stator leading edge are aligned azimuthally at the hub. The stagger angle of the vanes 
is ct s . The spacing between the blades is 27 rr / B, and between the vanes, 27T r/V . The local 
semi-chord of the vanes is b which can be a function of r. 
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CHAPTER 3 


MEAN ROTOR WAKE 

In this chapter, we provide information that will be needed for computing the acoustic 
response discussed later in Chapter 4. Our basic geometry will be a constant-area annular 
duct, through which a uniform flow is assumed. We shall describe a rotor wake model 
from which we derive an expression for the wake velocity normal to the stator vanes; from 
this expression we obtain the unsteady loading along the stator vanes. The loading will 
then be used in Chapter 4 to derive an expression for the complex amplitudes of the duct 
modes excited by the interaction of the rotor blade wakes with the stator vanes. These 
amplitudes will in turn give us the power propagating up and down the duct. 

3.1 Rotor Wake Model 

The pressure distribution on the stator vanes can be calculated by dividing the stator 
into a series of radial “strips” at constant radius and calculating the pressure on each 
“strip” as though it were a linear cascade of thin flat plates. In this approximation, the 
radial variation of the inflow to the vanes is ignored - - an application of “strip theory” 
to stator aerodynamics. We are concerned here only with the mean value (time average 
in rotor coordinates) of the rotor wakes, which, being time periodic in stator coordinates, 
generate sound at harmonics of the blade passage frequency Q.B /2tv. 

Consider a cylinder of radius r, centered on the duct’s axis of symmetry. The inter- 
section of this cylinder with the rotor and stator is depicted in Fig. 3.1. This surface is 
the one viewed by looking down the positive X 3 -axis in Fig. 2.1 toward the origin. In 
Fig. 3.1, two sets of coordinate axes are shown; axes (Xi, X 2 ) are fixed at the trailing edge 
of the rotor, while axes (xi, X 2 ) are fixed at the leading edge of the stator. The azimuthal 
coordinate, (r<j>), is in the opposite direction to the X 2 -axis as indicated in the figure. Rotor 
rotational speed, rO, is in this same direction. The (xi, ^-coordinate system is defined at 
each radius and places the leading-edge points of the stator vanes at the points (0, —vh), 
where v is any integer and h = 2tt r/V. As seen in the figure, u increases in the direction 
of increasing X 2 (decreasing </>). 
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The relation between the two sets of axes in Fig. 3.1 is 


X = x + D + Qrte 2 , (3.1) 

where <§2 is the unit vector in the azimuthal direction and D = D(r) is the vector distance 
from the trailing edge of a rotor blade to the leading edge of a stator vane. D is shown 
for the case where, at t = 0, the trailing edge of the rotor is aligned azimuthally with the 
leading edge of the stator at the hub (i.e., r = r H ). We will perform our initial analysis 
based on this configuration. However, D actually is defined only to within an additive 
vector r<foe 2 , where 4 > 0 represents an arbitrary angle of rotation of the rotor. We shall 
later indicate the effect of such a <f > 0 and assign it a particular value. 

In Fig. 3.1, W, which we will call the wake velocity downstream of the rotor, represents 
the entire air flow, including the mean flow, downstream of the rotor and in the rotor 
reference frame. The parameter a CL is the angle that the wake velocity W makes with 
the Xi -direct ion at the stator leading edge. 

So as to arrive at a reasonably simple model of the rotor wake, which will be valid in 
the vicinity of the stator vanes, it is convenient to introduce two plausible assumptions. 
They are 

1. no radial flow occurs, and 

2. pressure gradients and turbulent and viscous diffusion can be neglected 
over the chord of the stator. 

A result of the second assumption is that the wakes do not change during convection 
across the stator. The fluid velocity then has only two components, (Wi, W 2 ) in (Xi, X 2 ) 
coordinates, and the equations of motion of the fluid reduce to 

(W • V)W = 0, (3.2) 


V • (poW) = 0, 


(3.3) 



where />o is the nominal fluid density. A parallel flow of the form 

W = W(X ■ N) w, (3.4) 

where N and w are constant orthogonal unit vectors (see Fig. 3.1), is a solution to 
these equations, and is the form we will assume for the wake. The velocity is everywhere 
parallel to the unit vector w = (cosa CL ,sin o CL ), and is constant in magnitude on the 
lines X • N = constant. The variation of the magnitude across the wakes is determined by 
the as-yet-unspecified function W (X • N). 

In this report, our Fourier series are defined in terms of the pair 

OO 

m- £ ( 3 - 5 ) 

S= — GO 

X P /2 

f, = — f f(x)e- 2 ’’“ x/x 'dx, (3.6) 

Xp J 

— Xp/2 

where Xp is the period. Let H be the gap between rotor blades, measured azimuthally. 
Now, W(X-N) is periodic in the N direction, because all the wakes are assumed to be iden- 
tical. The period here is the normal distance between wake centerlines, Xp = H cos a CL . 
W{X ■ N) can therefore be written as a Fourier series, i.e., 

W{X ■ N) = jr e i [ sS ^ /(rcos “^ ) ], (3.7) 

S = — OO 

where we have substituted 2-jrr/B for H. The unit vector N has components (-sma CL , 
cos a CL ), so the Fourier series may also be cast in the following form, 

W = ^2 (3.8) 

OO 

which indicates that the rotor wake velocity is periodic in the azimuthal direction also (as 
well as the direction normal to the wakes). The azimuthal period is, of course, 2nr/B. 

In principle, the Fourier coefficients W s of the wake flow could be determined ex- 
perimentally by processing a sufficient number of flow measurements collected from a 
transducer mounted behind the rotor. If only the first few harmonics of the blade passage 
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frequency are of interest (it will transpire that the sth harmonic in the wake flow interacts 
with the stator to generate sound in the sth harmonic of the blade passage frequency), 
this is a feasible proposition. However, by assuming an explicit (and plausible) form for 
the wake velocity profile, it is possible instead to characterize the wake velocity profile 
semi-empirically by using experimental data to estimate the unknown parameters which 
appear in these forms. In practice, this data is taken at the stator leading edge; however, 
it must be referenced back to the rotor trailing edge to be appropriate for Eq. (3.7). In 
ri pfinin g the W,’s using Eq. (3.6), note that x is zero at the wake centerline for V072 .* 

There are presently three wake velocity profiles to choose from in V072: (i) Hyperbolic 
Secant, (ii) Gaussian, and (iii) Loaded Rotor 67. Three types of wake width and velocity 
correlations are available for use with these profiles: (i) Loaded Fan Wake Function and 
two types of (ii) Linear Rational Functions. These options will not be discussed here. For 
information regarding the wake profile and correlation treatment, the reader is referred 
to Refs. 16 and 17. In addition to wakes, the current version of V072 now also contains 
models for hub and tip vortex flows. These empirical vortex models, which were not 
part of the original code, were added later to provide a more complete representation of 
the flow field downstream of the fan. In their current form, however, the vortex models 
are intended, primarily, as a tool for providing crude estimates of the contribution of 
vortex flow to the rotor/stator interaction noise. As such, they should only be used in 
parametric assessments of hub and/or tip vortex noise contributions, not as design tools. 
When selected, these vortex models are combined with the wake profiles so that both are 
used in determining the W 3 ’s multiplied by an additional factor sin(a 5 + a CL ). This is 
the quantity that appears later in Eq. (3.18). These vortex models are discussed in 
Ref. 16. 

3.2 Stator Upwash Velocity 

We next deter min e the component of wake upwash velocity perpendicular to the sta- 
tor vanes, using the combined wake/vortex profile in rotor-fixed coordinates at the rotor 

* This is not the case for Ref. 16, where, with the notation here, x would be %p/2 at 
the wake centerline. This means that the W 3 ’s used in V072 differ from those defined in 
Ref. 16 by a factor of exp(is7r). 
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trailing edge, as given by Eq. (3.8). This upwash is specified in terms of its Fourier series 
harmonics by using a coordinate transformation to stator-fixed coordinates. 

We first calculate the wake velocity relative to the stator vanes, 

U r = W -flre 2 (3.9) 

and then find the component w of this velocity that is normal to the vanes, 

w = U r -h, (3.10) 

where h is the unit surface normal. Evaluating Eq. (3.10) using Eqs. (3.8), (3.9) and the 
coordinate transformation (3.1), we find that 


w 


°° r l 

= - Y, W * sin (a s +a CL )e-n sB(I?2 -' Dltantt ^ ) H 


(3.11) 


x e -i [ s - B ( X2 ~ Xltana cx,)/ r + sBnt ] 

where the W 3 ’s are the same as before, except that for s = 0, the Or part of Eq. (3.9) 
would be needed in evaluating Wo-* 

On any specific vane, say vane v, the components of x = (x\,X 2 ), can be found from 


x = zc + uh + 6c, (3.12) 

where c = (cos a s , — sin a s ) is a unit vector directed along the vane chords, and h = (0, h) 
is the vector separation between any two neighboring vanes. The coordinate z is equal to 
—b at the vane leading edges and to +b at the trailing edges. In terms of components, 
Eq. (3.12) can be written as 


X\ = (z + b) cosa^, (3.13) 

x 2 — ~ (z + 6) sin + 2i rur/V. (3-14) 

* Note, however, that the s = 0 case is not relevant to the work here since it represents 
the steady part of the vane loading. 
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Also, substitute for D\ and D 2 the relations 


— X SPAC^ 

^2 = Vrd “ 1 “ VsD 5 

which are easily obtained using Fig. 3.1. Then, Eq. (3.11) gives the result 


w = 



w e* \^ sZ “ ^ *] 


3 = — OO 


where 


w s = -W a sin(a <s + cn CL ) e 1 [ sB {Vrd + VsD I s'’^ tano cJ/ r ] 

x e* [ sBb ( sin a s + 008 “s tan a ci)/ r ] 


and 


(3.15) 

(3.16) 


(3.17) 


(3.18) 


cyjLJ . . v 

k a = — (smOg + cosa s tana Ci j. 
r 

In these equations, the parameters W a , a s , a CL , y RD , y SD , x SPAC , 
of the radius r. 


(3.19) 

and b are all functions 


Note that, had D in Eq. (3.1) included the additive vector r</> 0 e 2 , we would now 
multiply the right-hand side of Eq. (3.18) by a factor of exp(— isB<f>o). This vector, 
0o = —Xspac.h tan a crH /r, is used i 11 V072. Therefore, Eq. (3.18) becomes 


w 3 = -W s sin(a s + a CL ) e - i [ aB ^RD^ysD~ x sPAc t&na cL + x spAc,H tana ci„jf)/ r ] 


x e 


t [sBfc(sin a £ + cos at s tana ci )/r] 


(3.20) 


where the subscript H in Eq. (3.20) refers to the stator leading edge location at the hub. 
Thus the phase in the circumferential direction is defined such that the exponential part 
of Eq. (3.20) is unity at the stator leading edge location at the hub. 


3.3 Loading on Stator Vanes 

Because the vanes are impermeable, they must induce an additional velocity to cancel 
the normal inflow velocity at the vanes themselves. Thus vane pressure loading is generated 
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by the negative of w and is found numerically by solving an integral equation. This integral 
equation is obtained by modifying Eq. (A.21), which is derived in Appendix A. Eq. (A.21) 
can be written as 

+vbhi,X 2 + vbh 2 ) = e wa f K c {x i - y, x 2 ) y > ( 3 -21) 

-b ~ 1 

where w is the upwash evaluated at the point (xi + isbh \ , x 2 + ubh 2 ). The variables Xi, 
x 2 , hi and h 2 are defined in Fig. A.l, and u, po and b were defined previously. Parameter 
<j is the inter-blade phase angle, —2 irsB/V, and K c is the cascade kernel function, both 
of which are described in Appendix A. Note that the negative sign for a is a result of our 
convention for counting blades. Variable U r is the nominal fluid velocity in the rotor-fixed 
coordinate system (see Fig. A.l), and Apo is the pressure loading on vane v — 0, which 
has arbitrarily been selected to be the reference vane; Ap*, would be the loading on an 
arbitrary vane. From the discussion in Appendix A, it is easily seen that 

Ap„ = Ap 0 e^ = Apoe- 21 ™ 3 - 8 ^ 

gives the pressure loading on the i/th vane. 

If we restrict ourselves to vane v = 0 by talcing X\ — z and X 2 
where z is the local chordwise coordinate used before, we find that 

-(-0 (3.23, 

—6 

If w above is taken only to represent a single harmonic, then we use w s exp (ik 3 z) (for vane 
v = 0) from Eq. (3.17) in its place, and Eq. (3.23) becomes 

w 3 e ikaZ = J K c (z - y) y • (3-24) 

-b 

In Eq. (3.24), we have suppressed the 0 in the argument of K c and replaced Apo with 
Apo, 3 , where Apo, s is the sth Fourier series harmonic for Apo- Finally, we can write 


(3.22) 

= 0 in Eq. (3.21), 
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Eq. (3.24) as 


if we take 


+b 

e ikaZ = J K c (z — y) f a (r,y) y, 

-b 


(3.25) 


fs(r, z) 


poU r W a ' 


(3.26) 


In V 072 , we solve Eq. (3.25) for / s (r, z), the elemental stator vane loading function, to 
find the vane unsteady pressure loading. Given / 9 (r, z), the sth harmonic of the pressure 
on vane u, A p„ >3 , is given by 


= poU T w,f, (r, z)e~ !li "'- B ' v . (3.27) 

Relation (3.27), which will be used in the next chapter, follows easily from Eqs. (3.22) and 
(3.26). 

The parameters needed to evaluate K c , when solving Eq. (3.25), are listed below: 

/ 3 r : 

Reduced frequency: K = —sBVtb/ (/3^co), 

Inter- vane gap: h = ‘I'KrjV , 

Vane stagger angle: a s , 

T: -^^[1 + M T M r (r /ro) sin a s //3t]. 

These values appear in the discussion in Appendix A. In the above expressions, M r is the 
Mach number for the flow relative to the vanes, and M T is the rotational Mach number at 
the rotor blade tip. 
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CHAPTER 4 


DUCT ACOUSTICS 

To calculate the duct mode amplitudes excited by the interaction of the stator vanes 
with the mean rotor wake, we first have to introduce the notion of normal modes. We will 
then use the loading on the stator vanes from the wakes, in combination with a Green’s 
function integral, to obtain the modal amplitudes. From these modal amplitudes we derive 
expressions for sound power in the duct. 

4.1 Normal Modes in an Annular Duct 

The normal modes of an annular duct with hard walls are the set of solutions to the 
two-dimensional Helmholtz equation having the special form 

*(r»=/(r)e 4m ^, (4.1) 

where m is any integer, positive, negative, or zero. The Helmholtz equation is given by 


V 2 # + k 2 '& = 0, 


where V 2 is the two-dimensional Laplace operator, in polar coordinates, 

2 d 2 Id Id 2 

V = + + (4 ' 3) 

In this equation, k is an undefined constant. As is shown below, the equation has non- 
trivial solutions which fit the appropriate boundary conditions on the walls of the annular 
duct only for certain specific values of k. If Eq. (4.1) is substituted into Eq. (4.2), we 
obtain 

+ I ^) / = 0. (4.4) 

dr z r dr r l 

Substituting u — k r for r reduces Eq. (4.4) to Bessel’s equation of order m: 


fl + I 

du 2 u du 


1 df (, 

udu V 


/ = 0 . 
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The solutions axe mth order Bessel functions of the first and second kind, i.e., f(u) = 
^> m (/cr), where 


i>m(Kr) = AJ m (Kr) + BY m {n r). (4.6) 

If we specify that 4»(r, <f>) — ip m (n r) exp is either the pressure or the velocity 

potential, then the radial derivatives of rp m must vanish at the inner and outer walls of the 
duct, so if A and B are to be not both zero, 


J'm( Kr D) 


Y m( Kr H ) 

Y L{^ r D ) 


= 0 . 


(4.7) 


This transcendental equation has a countably infinite number of roots k for every integer 
m; if we denote these roots by K mn , n = 1, 2, 3, • • • , and arrange them in order of increasing 
magnitude, then the functions V>m(«mn r) have (n - 1) zeros in the interval r H <r <r D . 


These functions are also orthogonal with respect to the weight function r over the 
same interval; that is, 


' D 

I 




(4.8) 


unless i — n. Using this fact, plus the fact that the functions exp m = any integer, 

are orthogonal on the interval 0 < <j> < 2n, it is easy to show that the normal modes 
ipm( K mn r ) ® X P (twi^) are orthogonal over the cross-section of the duct. Specifically, 


27T 


// 


{^m(«r 


r ) e im<t> } x {^(/tfc/r) e tk<t, }rd<j>dr 


0, if k ^ m or i n; 

r D 

J ipm( K rn n r) rdr, if k = m and £ = n. 


(4.9) 


It is convenient to adjust the constants A and B in Eq. (4.6) (which are determined by 
Eq. (4.7) only to within an arbitrary multiplicative factor) so that 
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(4.10) 



The orthogonality relation, i.e., Eq. (4.9), then becomes 


r D 2tt 

J J {lpm(K 


mn 


r) e lTn4> } x {ipk(Kker) e lk<t> } rd<f>dr = w(r 2 - 


r H 0 

where 6 m k and 6 n g are Kronecker deltas; S m jc is equal to unity for m 
for m ^ k. 


rZ)6mkt W, (4-11) 

= k and equal to zero 


A special FORTRAN program has been written which calculates K mn and the constants 
A and B in Eq. (4.6) subject to the normalization given in Eq. (4.11); for details see 
Appendix B. 

The significance of the normal modes is that they can be used to represent pressure 
patterns which propagate within the duct without change in form, and that any acoustic 
field within the duct, however generated, can be represented as a suitable combination of 
these patterns. 

The wave equation in a duct containing a fluid moving at a uniform axial velocity U 


2 d 2 p 1 / d d \ 2 

Vp+ M~^(9t +£ W P = °’ 


(4.12) 


where p is the acoustic pressure, and V 2 is the two-dimensional Laplace operator, defined 
by expression (4.3). 

Now assume that p has the form 


Ipm(Kmnr) S”*' ~ ^~ sB * , (4-13) 

where we have assumed for frequencies the blade-passing values, sBfl, introduced in the 
last chapter. The pressure pattern above has m diametral nodes (like the spokes of a 

* Between the rotor and stator an appreciable circumferential velocity exists as well. 
T his velocity is not accounted for in Eq. (4.12). 
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wheel), and (n - 1) concentric circular nodes. If we substitute expression (4.13) into Eq. 
(4.12), and recall from Eq. (4.2) that 

V 2 [*,(«„ r) e im *] = -K? mn lp m (K mn r) e im +, (4.14) 

the following relation is obtained: 

-L-7’ + (f^) 2 =°. (4-15) 

Solving for 7 in terms of Kmn and sBft, we set 7 = 7 mn s and find that 

7mnj = ~jj 2 y “ (4.16) 

where the plus sign is used for upstream propagating waves and the minus sign for down- 
stream ones. Also, 

k m „, = y(^) 2 -/? 2 <„. (4-17) 

In the above expressions, M is the axial flow Mach number, U/cq, and (5 = y/l - 

Eq. (4.16) shows that 7 mns is real whenever sBQ is real and \sBQ./cq\ > /3 Kmn ■ Under 
these conditions, the pressure pattern propagates unchanged along spiral paths normal to 
the lines 

m<f> - 7 mna = constant, (4.18) 

producing a rotating pattern. If sBfl is real but |sBfI/co| < ■> then 7 mns is complex 

and the pattern grows or decays exponentially along the duct depending on which sign in 
Eq. (4.16) is selected. 

Modal pressure patterns such as 


lpm{Kmnr) ^s*! ~ tit) (4.19) 

can be superimposed to form a general description of the acoustic pressure field in a 
turbomachine. Thus, if we multiply the pattern above by an arbitrary coefficient, say 
Vmnsi and sum over m, n and s, we obtain 
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(4.20) 


oo oo oo 

P (x, t)= y, EE iM«m» r ) ei{m * ~ ymnaxi _ sB n t} 

s — — oc m=— oo n=l 

where now a; is the point (r, Xi, <j>) in the duct. Note that in general p(x,t) is composed 
of both propagating and non-propagating modal pressure patterns. The Fourier inverse of 
this equation gives the sound pressure at the sth harmonic of blade passage frequency: 

OO OO 

P.( S)= 52 52 (4.21) 

m= — oo n= 1 

It should be noted that, following standard practice, when output is printed in the V072 
code, the radial mode index is adjusted to begin at 0 instead of 1. 

4.2 Duct Acoustic Modes 

The acoustic modes of the inlet or exhaust ducts were just derived above. In this 
section, we shall derive equations that relate these modes to fluctuating loads on the 
stator vanes. The end result will be an expression for the complex pressure amplitudes 
appearing in Eq. (4.21). If we assume a harmonic time dependence for the excitation, say 
f s (x) exp(— isBCtt), which is the case for the wake, then pressure at the field point x will 
also be harmonic with the same frequency, and will have the form p s (x ) exp(-isBQt). 
This is due to the linearity of the governing equations. 

It can easily be shown* that the pressure fluctuation p s (x) within the duct is given, 
in terms of the force/unit area f s (x ) exerted by the vanes on the fluid, by the Green’s 
function integral 


?.(*)- J {V f G(f,j 7)} ■ f. (S) dS(ff), (4.22) 

S(y) 

where S(y) represents both surfaces of the vanes and Vy is the gradient operator with 
respect to the vector y. In this equation, G(x,y) is the space-only dependent Green’s 

* One way to see this is to start with Eq. (4.13), the time dependent Green’s function 
integral for density fluctuation p(x, t), in Ref. 19. Multiply through by Cq to switch to 
pressure fluctuation p(x,t). Then proceed much as is done in section 4.3.2 (“Application 
to Pure Tones”) of Ref. 19 in obtaining Eq. (4.25) there. 
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function for a hard- walled annular duct. It is the same function as G w (y \ x) given by 
Eq. (1.C.12) in Ref.19. If we introduce a cylindrical coordinate system such that the field 
point x is (r,Xi,<(>), while the integration point y is {r' ,y\, </>'), then G can be written as 
follows 


^ 3L ' l Pm.{. K mn r ) *i)] (4 23) 

2 i , ^mni r 


m=—oc n= 1 


where T = 7r (r^ — r£). 

The coordinate system used is pictured in Fig. 4.1. This figure is obtained by con- 
sidering an imaginary cylinder, of radius r, centered at the duct axis of symmetry. The 
intersection of this surface with the stator is what is depicted. The origin is located where 
the leading edge of the stator intersects the hub. In contrast to the situation in Fig. 3.1, 
the origin is always here rather than moving with the leading edge as the radius changes. 
In expression (4.16), which defines 7 mR „ the plus sign is used for upstream field points 
(xi < j/i) and the minus sign is used for downstream field points. 


If the fluid is assumed to be inviscid, the force f a will be normal to the surface, 


fs=P,n, (4.24) 

where n is the outward surface normal and p 3 is the local pressure. Thus 

p 3 (x)= j n(y) {VyG(x,y)}p a (y) dS{y). (4.25) 

s{y) 

The integration is to be carried out over both faces of each vane. On each vane, we can 
divide the integral in the equation above into two parts, one over the forward or upstream 
face of the vane, and the other over the downstream face. Denoting the upstream face by 
superscript (+) and the downstream face by superscript (— ), we have, because the vanes 
are very thin, 


P a n = (p a - pt) h, ( 4 - 26 ) 

where h is the unit normal vector erected on the mean surface of the vane as shown in 
Fig. 4.1. Define A p a as A p 3 = p~ — pf , which is positive when the pressure is greater on 
the lower face of the vane. Then 
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(4.27) 


P,(x)=-J ft®- {VyG(l,t/)} Ap,® dS(y), 

where now the integration is over the vane mid-surfaces S M (as opposed to the two vane 
faces). 

Using expression (4.23), Eq. (4.27) can be written as 


oo oo 

p,® = E E 

m=-o° n=l ( 4 .28) 

x— i / ft® -V 5 [€«-"*' + ->•-*>")] Ap,®dS®. 

Jll rCrnns J L 

Note that ^m^mn^) exp[i(m0 - 7 mns Xi)] is a rotating pressure pattern of the type 
discussed previously.* Eq. (4.28) is thus a normal mode expansion of the sth harmonic of 
the acoustic pressure within the duct. Comparing Eqs. (4.21) and (4.28), we see that 


Pmns 


2irfc r 


J Ipm(Kmnr') fl{y) • Vy [< 


,l(— m<t> 1 -f 7mnS Vl) 


A p 3 {y) dS(y), (4.29) 


where the p m ns are the complex modal amplitudes introduced in Eq. (4.20). Eq. (4.29) is 
the starting point for evaluating these coefficients. 


In Eq. (4.29), the integration is to be carried out over all V stator vanes. As before, to 
reduce the region of integration to one vane, we arbitrarily select one vane as a reference 
vane. We assign this vane the number 0, and the remaining vanes, in the direction of 
decreasing <f>* , the numbers 0 through {V — 1). Using the notation in Chapter 3, let the 
pressure on the vth vane be Ap^yo), where y 0 is the point (r',y 1? ^ / ) on the v = 0 vane. 
Then the corresponding point on vane v is (r / ,yi,0 / — 27 tis/V). The total contribution of 
all V vanes is thus 




Pmns = 7PFZ / ^rn{^mnr') fl{y Q ) ■ Vjf 

2,11 K m ns J 

So 

x Ap„,,(!/o) e 2 "”“' /V } ASo(po), 


(4.30) 


* Except for the factor exp(— isBSlt). 
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where we have gone from a notation where S M represented the middle surfaces of all the 
vanes to one where So represents the middle surface of the reference vane. 

Referring to the geometry in Fig. 4.1 and noting that the normal vector on the 
reference vane is n = (sin a s , cos a s ), we have that 


n • Vjj- e i ^~ Tn4> ' + ’ rrnnsyi ' ) = i (^7 cosa s + 7 mn 3 sina s ) m<t> +f mns y i '>. 


(4.31) 


Note that because <f> and y 2 are in opposite directions, there is no minus (-) sign preceding 
the m/r' term. To facilitate the integration over the vane surface use the intrinsic chordwise 
coordinate z 7 , which as before varies from —b at the leading edge to +b at the trailing edge. 
From Fig. 4.1 we see that 


yi = - x SD + b cos a s + z 7 cos a s , (4.32) 

<t> = (- y SD + bsina s + z 7 sina s )/r 7 . (4.33) 

Using r' and z' as integration variables, and applying Eqs. (4.31) - (4.33), Eq. (4.30) 
becomes 


D 

Pmns ~ TypT / r ) COS QC S + 7 mns s ' m a s J 

ZL Kjxxns J 'T 

r H 

r V_1 

xc i(-^n«.* SI ,+ m, SJ ,/r') n Ap^r^z') e 2nimu/v ] 

- b 1 *'=° 

x gi [bmns cos a s - m sin q s /r')(z'+t)] dz' dr' . 


(4.34) 


The sth harmonic of the pressure loading on vane u is given by Eq. (3.27). If we substitute 
this into Eq. (4.34), we have 
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Pmns — 


P 0 


2rfc„ 


'D 

J ipm ( K mn r') U t w 3 cos oc s + 7 mn3 sina 5 ) 


x e 


0 2 > nii/{m—sB)/V 


} 


V-l 

'TmnS X SD m ^SD / r ) | ^ ^ £" 

*/=0 
b 

X / T,[r\Z ) jli-rmnscosas-msmacs/r') (z'+6)] ^ rfr ' 
-fe 


(4.35) 


o 

J fa{ r 'i z> ) 


The sum over vane number in Eq. (4.35) can be evaluated explicitly: 

V, if m — sB = — gV'; 


v-i 


^ ^27ru/(m— sB)/V _ 

t/=0 

Here q is any integer. The final result is 


-r 

U 


e 

=0 1 0, otherwise. 


Pmns — n-pil f 1pm{ K mn r ) U r W 3 (r) ( COSQ! s + 7 mns sin Q: s ) 

Zl Kmns ' 

i, (4.36) 

x e i(— + f fJr, Z ) dzdr, 

-b 

where now the primes have been dropped from z and r, and the index m is given by 

m = sB — qV. (4-37) 

Eq. (4.36) is our final result. As is evident from Eq. (4.36) above, p m ns is propor- 
tional to w 3 , the sth-order Fourier coefficient of the mean rotor wake. The computer code 
computes these coefficients, as mentioned previously, using semi-empirical information. 
However, the code can easily be modified to accept the Fourier coefficients of the wake 
velocity as inputs, and compute the duct modes directly from them. Additionally, new 
semi-empirical wake modules could be added to V072 to supplement those already there. 

As a final remark, let us note that once the pressure amplitudes, p mnJ , are specified 
by virtue of Eq. (4.36), the pressure at any point within the duct can also be determined, 
if desired, from the summation 


- 30 - 



oo oo oo 

p(x,t)= Y, E £ (4.38) 

q= — oo n=l 5= — oo 

where m = sB — qV. 


4.3 Sound Power 

The flux of sound power in the duct, is found by substituting Eq. (4.36) into Eq. 
(C.20) of Appendix C, which is written as 


Power = 


P oU 


OO OO OO 


EEE 

q— — OO 71=1 8— — oo 


(4.39) 


Here G m na is defined as 


Gmns — 


qFikf 2 /? 4 ( sBn/U)k mr 


(4.40) 


(sBCI/cq ± Mk mn3 ) 2 

The upper set of signs apply upstream of the stator and the lower set apply down- 
stream. The details describing how power is determined are presented in Appendix C. 
In Eq. (4.39), the index m is given by Eq. (4.37). Finally, if the sound power at one spe- 
cific harmonic of the blade passage frequency is desired, the summation over s is deleted, 
and s is set equal to the desired harmonic number. 


Note that for a propagating mode (m,n) at harmonic s, another mode (— m, n) of 
complex conjugate amplitude propagates with frequency — s. Because the amplitudes are 
complex conjugates, V072 only calculates positive harmonics. The total sound power for 
the mode is then evaluated from just the one by applying a factor of 2. 
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CHAPTER 5 


NON-DIMENSIONAL RESULTS 

The V072 program uses non-dimensional input geometry and performance parameters. 
It uses these quantities to compute dimensionless versions of pressure amplitudes and power 
and then converts the results to dimensional values. The final output of the program is the 
magnitude and phase of each propagating mode excited by the interaction of the stator 
vanes with the wake of the rotor. The sound power flux per mode is also computed and, 
by summing over all propagating modes, the total sound power flux is obtained. Note that 
propagating modes are those discussed in Section 4.1 for which sBQ / Co > ft K mn ■ Modes for 
which sBQ/co < pKmn are said to be cut off and are not included in the output. Note also 
that a multi-vaned stator, as seen by Eq. (4.37), excites only a subset of the propagating 
modes at any given frequency; specifically, only modes whose number of diametral nodes, 
m, is related to the number of rotor blades and stator vanes by the equation m = sB — qV , 
where s and q are arbitrary integers. The program takes this selection mechanism into 
account in choosing which mode amplitudes to compute. 

In perfor min g its calculations, V072 is divided into two parts. It first carries out the 
rotor wake calculations and then computes the acoustic results. To do this, it needs both 
geometric and performance input. Note that performance data would normally come from 
a steady aerodynamic prediction code of axisymmetric or two-dimensional nature. (See 
Ref. 3 for more discussion of this.) 

Non-dimensionalization in V072 is with respect to a reference pressure ^pqU 2 , a ref- 
erence power gpoCoM 4 rJ, lengths r D and b T , and speeds cq, U and U r . Here po is the 
nominal flow density, U the mean axial fluid velocity, Co the speed of sound in the duct, 
M = U / Co, r D the duct outer radius, and b T the vane semi-chord at the tip. 

5.1 Pressure Amplitude 

The dimensionless expression for the complex pressure amplitudes is easily derived 
starting from Eq. (4.36). After straightforward effort, it takes the form 
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Pmns 


a r V 


\poU 2 7r(l - a%)k r 


f 1 fU r \ 2 . . ,m _ . x 

J \~u) <bs ^ m ( Xmn x > cos a 5 + 7mn» Sin ) 

X C i[(-W s +my s /i)<r c ] 6 C mn3 (x) (fo. 


(5.1) 


In the above, a c = 2 b T /r D \ a r = r^/r^ where r H is the hub radius; 7 mn9 = 7 mn3 r D , where 
jmns, the axial wave number, is the same as in Eq. (4.16); and k mns = k mns r D where 
kmns is the same as in Eq. (4.17). Note that 7 mn3 and k mns can be written, explicitly, as 


7 mnj ^2 ± imi 


(5.2) 


and 


fc mn3 = (5-3) 

where M r is the rotor rotational Mach number at the blade tip and j3 = \/l — . Further, 

in Eq. (5.1), X mn = Km n r D where K mn is the eigenvalue for mode (m,n); x = r/r D ; 

x s = x sd /^t > Vs V sd /^t > ^nd ^ 6/26 t . 

Additionally, C rrm3 (x) is the same chordwise integral of the elemental blade loading 
function f s * multiplied by an exponential function, as seen in Eq. (4.36), only here non- 
dimensionalized using the chordwise parameter z/b. It is given by the expression 


Q _ j j 3 e i [ff c i(7mns cos -m sin a s /x) (l+z/6)] 


(5.4) 


Because / 3 has a square root singularity at the leading edge (z/b = — 1), it is necessary to 
change the integration variable from z/b to ip, where z/b = cos ip. We then have 


Cmnsix) = fs Sin if e [-cKWs cos -m sin A) (1+cos ^)] ^ ( 5 .5) 

JO 

This integral is easy to compute numerically using Simpson’s rule, because the product 
f s sin ip is finite at the leading edge (ip = 7 r). 

* Note that the elemental loading function / 3 becomes a different function as its inde- 
pendent variable changes from step to step. However, for simplicity of notation, we still 
denote it as / 3 . 
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The parameter w a , used in Eq. (5.1), is defined as w 3 /U r , where w a is the quantity 
given by Eq. (3.20) and U r is given by Eq. (3.9). It can be written as 

W = W e~ iSB l aC ( y * +y S )/X ~ X SOR t * Ua CL+ X SOR,H tana CL,H\ 

9 9 r - , (5.6) 

^ isB L<r c 6 (sin <x s H-cosa^ tana CL )/xJ 

where 


tis = ~ Sin(a s + a CL ). (5.7) 

Also, y R Vrd/’^T % SOR SPAC 

The spanwise integration in Eq. (5.1) requires special treatment because w s , as seen 
in Eq. (5.6), contains a phase angle, x SOR tan a CL , which varies rapidly over the span 
when the separation between rotor and stator is large or when the relative flow angle at 
the stator leading edge, a CL , is a strong function of the radius, r. To handle this rapidly 
varying phase, an adaptation of Filon’s rule of integration has been used for the spanwise 
integration (see Appendix D). Note also that for spanwise integration, V072 requires values 
of geometric and performance quantities at radial locations determined internally by the 
program. To obtain these parameters at the necessary locations, V072 linearly interpolates 
the pertinent input data. 

To obtain f a (y/b), which is needed in Eq. (5.4) to evaluate C rnna (x), one solves the 
dimensionless version of Eq. (3.25) which is easily seen to take the form 

e ik.z/b _ f 1 (5.8) 

Jy/b =- 1 b 0 o' 

whose kernel function K c ((z — y)/b), as mentioned previously, is derived in Appendix A. 
The quantity k a is the non-dimensional version of the vane chordwise wavenumber k s (see 
expression (3.19)). It is given by 

- sBa c b . . . . . 

k a = k a b = ( Smc *s +cosa s tan<* CL ). (5.9) 

Finally, there are issues as regards a singularity at the leading edge. For this reason 
Eq. (5.8) must be changed to the form 
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(5.10) 


K c (cos6 — cos ip) F a (ip) dip, 

using the same chordwise type transformations as used previously, i.e., zjb — cos 6 and 
y/b = cos ip. In Eq. (5.10) we then have F a (ip) = f 3 (cos ip) sin ip. The singularity issues 
along with a solution of Eq. (5.10) using the method of collocation (Refs. 20 and 21) are 
discussed in Appendix E. 

Note that, in the code, magnitudes and phases of the pressure mode amplitudes are 
printed out. The phases are given for positive values of the index s. The magnitudes are 
r.m.s. values obtained by time-averaging the square of the sum of the -I -s and — s waves, 
{pmns 6xp( isBClt) + P(_ m )n(-s) exp(isBfit)} 2 . Since p(_ m ) ra (_ a ) = Pmns) the result 
of this calculation is y/2 | p mna |. Before output, the magnitudes are converted back to 
dimensional form through multiplication by ^ poll' 2 . The final result is then presented in 
dB units relative to a reference pressure of 0.0002 x 1.4504 x 10 — 4 psi (i.e., 2 x 10 5 Pa). 


„ik s cos 9 


= / 

JO 


5.2 Sound Power 

Given the modal amplitudes Pmns/ {\poU 2 ), the sound power is obtained from Eq. 
(4.39). Non-dimensionalizing by ^poC^M^r^, we have 

Power _ 27T (r£—r£)Gmns 

WlM^l ~ ^ Mr* 

where G mna is the quantity given by Eq. (4.40), m = sB — qV , and sound power in Eq. 

(5.11) is for one specific harmonic, s, of the blade passage frequency. In obtaining Eq. 

(5.11) , we have used the fact that \pqCqM a t^ is equivalent to \poU z Mr^. Eq. (5.11) can 
finally be re-written as 


Pmns 


ion !/ 2 


(5.11) 


where 


Power _ ^ 

ipocBMVJ = 


V 2jt (1 - aj) G mn , 

n=l 


Pmns 


ipoU 2 


(5.12) 


G m „, = T bum sJMtP (5.13) 

(sBM t ± Mkmna) 

Again, as for expression (4.40), the upper set of signs apply to power for waves propagating 
upstream and the lower set apply downstream. 
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V072 gives values of sound power for three different cases: (1) for the sum of all modes 
with a given s; (2) for the sum of all modes with a given s and m; and (3) for individual 
modes with given s, ra, and n. Case (1) is treated using Eq. (5.12); case (2) using Eq. 
(5.12) with the summation over q deleted and q set to give the value of m desired; and 
case (3) by taking the individual term in Eq. (5.12) corresponding to m, ra, s. For all 
three cases, as for the mode amplitudes, the actual power output comes from summing the 
values for both +s and -s. Additionally, the power is converted back to dimensional form, 
through multiplication by ^po<^M 4 r^, and then given in dB units relative to a reference 
power of 10 -12 /1.3558 ft-lbs/sec (10 -12 watts). 
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CHAPTER 6 


CONCLUDING REMARKS 

This report is an updated version of the analytical documentation for the V072 Rotor 
Wake /Stator Interaction Code. It provides updated geometry and revised equations and 
has eliminated the material related to turbulent acoustic sources, because these sources 
are not active in the code. Additionally, equations are now developed in terms of integer 
multiples of blade passage frequency, rather than the more general frequency u> used when 
turbulence was covered. 

In this report, equations have been derived for the amplitudes and power of the prop- 
agating duct modes excited by a turbofan (rotor/stator stage) operating at subsonic tip 
speed within an infinite hard-walled annular duct. The equations used in the code are 
non-dimensional and the output is tone noise generated by the mean velocity defect wakes 
of the rotor blades impinging on the stator vanes. Output is provided at blade passage 
frequencies for the propagating (cut on) modes. 
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APPENDIX A 


KERNEL FUNCTION FOR A LINEAR CASCADE 
IN SUBSONIC FLOW 

The kernel function for a cascade of thin airfoils in oscillating subsonic flow has been 
derived by several investigators, using a variety of methods (Refs. 5-9, 22). Of the vari- 
ous approaches employed to date, the method of Fourier transforms is perhaps the most 
straightforward. This procedure was used in Ref. 9, for example, but the inversion of 
the Fourier transform was accomplished numerically. More recently, Goldstein (Ref. 19) 
pointed out that the Fourier transform of the kernel function contains no branch points, 
so the inversion can be accomplished quite easily by using the Cauchy residue theorem. 
Goldstein outlined the procedure to be followed, but did not actually carry out the calcu- 
lation of the kernel function. The purpose of this appendix is to set forth the details of 
the inversion, and to record the end result. For the convenience of the reader, as well as 
to document the notation used, a brief derivation of the transform of the kernel function 
is also presented. 

The cascade geometry is shown in Fig. A.l. The airfoil semi-chord is b, and the gap 
between the neighboring airfoils is bh, with components bh x , projected along the chord, 
and bh 2 , normal to it. The airfoils are shown as having no camber, because the ultimate 
objective is to calculate the pressure field scattered by the cascade when it is subjected 
to vorticity convected with the mean flow. However, other situations, such as a cascade 
of oscillating airfoils, can be handled as well. For the purpose of calculating the kernel 
function, we need only suppose that a known chordwise pressure distribution exists on 
each airfoil, and calculate the resulting velocity field. 

The first step is to calculate the upwash generated by a single airfoil. Let the pressure 
be given by the real part of pexp(— isBSlt), and the corresponding velocity field be the 
real part of (w x , w 2 ) exp(-isBClt). Then, p satisfies the convected wave equation, and 
the velocity field is related to the gradient of p through the equations of motion.* These 
equations are 

* p and (w x , W 2 ) now represent the sth harmonics of pressure and velocity field. 
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(A.1) 


d 2 p 


d 2 P 

dy\ 


\(-isBQ + U rl ?-) 2 p = 0, 
c o &yi 


(-isBCl + U r -^—)wj + —Tr~ = 0; 3 = I) 2 - (A.2) 

v dyi Po dyj 

Consistent with the definition of Fourier series in Chapter 3, the following definition of the 
Fourier transform pair is chosen: 


+oo 

p= h I p eiayi<kx > (A - 3) 

“OO 

+oo 

p= j p e~ iayi dyi. (A.4) 

— oo 


The Fourier transforms of Eq. (A.l) and the second of Eqs. (A.2), i.e., for j — 2 are 


d 2 p 

dyl 


- X 2 p = 0, 


(A. 5) 


where 


i (— sBQ + a.U r ) W2 + 


1 dp 
po dy 2 


= 0 , 


(A-6) 


A = [a 2 - M 2 (a - sB£l/U r ) 2 ] h . ( A.7) 

The pressure p generated by an isolated airfoil is bounded at infinity, and, because the 
airfoils are assumed to have zero thickness, is antisymmetric in y^- A solution of Eq. (A.5) 
satisfying these conditions is 


p = - e A|ya| sgn(j/2), (A.8) 

provided that A is defined so that its real part is non-negative on the path of integration 
used to invert the Fourier transform. In Eq. (A.8), A p(a) is the Fourier transform of 
the chordwise pressure distribution, A p, on the airfoil. (Ap is positive if the pressure is 
greatest on the lower face of the airfoil.) By eliminating p between Eqs. (A.6) and (A.8), 
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we obtain the Fourier transform of the upwash generated by a single airfoil. The upwash 


itself is 


+ 0C 



— OO 


A 

2 i (a — sBSl/U r ) 


e *atyi— A| j/ 2 | d a 

POU? **' 


(A.9) 


The upwash generated by a cascade of airfoils, located at the points (j/i, 3 / 2 ) = m{bh\,bh 2 ), 
m = 0, ±1, ±2, . . is obtained by summing the contributions of the individual airfoils: 


» +OO 

W2 , v _ 1 f ^ Apm(a) 

u'r V2 ~~ 1 ^ Tn t^ 00 J 2i(a - sBCl/Ur) poU? (A.10) 

x e ia(y 1 -mbh 1 )-\\y2—mbh 2 \ ^ 

This infin ite series can be summed analytically if the transformed pressure distribu- 
tions on successive airfoils in the cascade [i.e., Ap m (a)] are related by a constant increment 
in phase angle. That is, for any integer m, 


A p m (a) = Ap 0 (a)e tm<r , 

where Apo(o:) is the transform of the pressure on the “zeroth” 
a (called the inter-blade phase angle) is a constant given by* 

27rsB 


(A.ll) 

(i.e., reference) airfoil, and 
(A.12) 


To calculate the upwash near any selected vane, say vane u, introduce 


x\=yi~ vbh\, (A. 13) 

5C2 = S/2 — vbhv, (A.14) 

and 

n = m — v. (A. 15) 

* This value can be derived by proceeding as in Ref. 23, page 38, only using our 
definitions of geometry and parameters. 


-41 - 



Then 


H)2 

— (x i + ubhi,X2 + vbhz) 

Ur 


+oo 

f Ap 0 (a) 
2t r J p 0 U? 

— OO 


\S(oc)e' ax ' 

2 i(a — sBSl/Ur) a* 


(A.16) 


where the infinite series 


can be summed as 


for 


«(“)= E ^(in<r— \\x2—nbhi\—ianbh2) 


S(a) 


1 ^^2) 

2 sinh(^A+) 


e (lA_ + Ax 2 )' 

sinh(^A-) 


(A.17) 


(A.18) 


0 < X 2 < 


(A.19) 


where 


A± = ±\bh ,2 + i(cr — abh\). 


(A.20) 


Now, use the convolution theorem to calculate the inverse of Eq. (A.16) 


yr{x i + vbhi,X2 + vbh2) = e* 

u r 


+b 

J K c (x i-(,x 2 ) 


Apo(C) <K 

PoU? b ’ 


-6 


where AT C is the desired cascade kernel function given by 


(A.21) 


+ OO 

. , b f A Se taXl 

^1,X 2 ) = — - J 2i(a - sBQ/Ur) da. 

— OO 

The factor b is included to make K c dimensionless. 

It is convenient to introduce dimensionless variables, as follows: 


(A. 22) 
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(A. 23) 


K = -sBQb/itf co), X! = x x /b, 
a = ab — M r K, x<i = x^jb. 

where j3 r was defined at the end of Chapter 3. The shifted transform variable a completes 
the squares in A: 

A = ^Vd2_K2. (A. 24) 

b 

Then 


where now 


K c = — 


0 r e 


iKM r x i 
2ni 


+ 0O 

/ 


'ySe 


,iaxi 


2 (a + K/M r ) 


doc, 


(A. 25) 


7 = (a 2 -if 2 ) 5 , 


(A.26) 




,(^+-^ 7 * 2 ) +/? r 7X 2 ) 


sinh(|A+) sinh(^ A_) 


(A.27) 


Notice that 7 S' is an even function of 7 , so that even though the integrand contains 
the variable 


7 = \fd? — A 2 , (A. 28) 

there can be no branch points at a = ±K. Thus, if we apply the residue theorem to 
evaluate the kernel function by closing the path of integration on a large arc in the upper 
or lower ci-plane, no residual integrals around branch cuts appear; the kernel function is 
simply the sum of the residues in the upper or lower half plane. To insure that the integrals 
on the arcs vanish as their radii are allowed to become infinitely large, the integrand in 
Eq. (A.25) must be modified. First, note that 


o 1 g!g 
(A- 7)2 9xl 


(A. 29) 
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so K c can be written as follows 


K c = 


giKM r x i 

2nif3 r 


_ +°° 

d 2 f Se totXl 

W 2 J 2 7 (o + K/M r ) da ' 

— OO 


(A. 30) 


The integral above can be evaluated via the residue theorem by using the contours 
shown in Fig. A.2, and the differentiation with respect to £2 carried out afterwards. If 
X\ > 0, the integral around the semicircle in the upper half plane vanishes when the radius 
of the contour goes to infinity. Thus, 



iKM r x 1 


f sum of residues in 
dx\ \ upper half plane 


(A.31) 


On the other hand, when X\ < 0 the integral around the contour in the lower half plane 


vanishes, so 


J_ JKM r *i d 2 f sum of residues in 1 
" j3 r dx\ \ lower half plane J 


(A.32) 


The integrand in Eq. (A.30) has poles at a — —K/M r , and at points where 
sinh(^A±) = 0. The latter set of points are located at ot± where 


A± = 2mri ( n = any integer) (A.33) 

[cf. Eq. (A. 27)]. When A± is real, some of these poles lie directly on the real axis. To 
arrive at the correct expression for K c , it is necessary to invoke the causality condition by 
stipulating that K has a small negative imaginary part. Once the kernel function has been 
evaluated, we can let Im(AT) — > 0. The effect of this procedure is to eliminate the possibility 
that acoustic waves not generated by the cascade itself are inadvertently included in the 
solution. 

If Im(A') < 0, the pole at & = —K/M r is clearly in the upper half plane. The residue 
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is 


r k 


M r 
2 PtK 


sinh (ffKh 2 /M r ) 

cosh (fiKh/Mr) - cos(r + Khi/M r ) 


—iKxx/Mr 

C 5 


where 


(A.34) 


r = a — abhi T7~- (A. 35) 

M r 

Using Eq. (A. 12) for cr, Eq. (A.23c) to obtain a when a = —K/M r , and Eq. (A.23a) for 
K, we can rewrite Eq. (A. 35) as 


1 + M T M r (r/r D )sma s /Pt , 
which is the expression for T used in the code. 

The roots ot^ of Eq. (A.33) are given by 


r = 


2nsB 

V 


(A.36) 


where 


a. 


T n hl Prh<2 

dP d 


[K 2 - (r n /d) 2 ]*, 


(A. 37) 


r n = r - 2n7r, (A.38) 

d = (A. 39) 

If the square root in Eq. (A.37) is defined as follows, 

[Re i6 ] * = \fRe i0/2 , 0 < 0 < 2tt, (A.40) 

where \/~R is the positive square root of R, then solutions having the plus (minus) sign 
in Eq. (A.37) are located in the upper (lower) half plane. In either case, the residue is 


Rn ( a n ) — 


j3 r h2 

2d? 


[(at -Tnhx/d?) (at + K/M r ) 


i 


(A.41) 


Using Eqs. (A.31) and (A. 32), and carrying out the indicated differentiations, we 
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obtain the final result 


K =&- 
C 2 Mr 


smh(l3^Kh2/M r ) 


cosh.((3? Khq/Mr) — cos(r + Khi/M r ) 


,-iPlKxi/M r 


(A.42) 


2d? 


+ oo 

E 


K 2 - * 2 ) 


L (at-Tnhx/d?) (< oct + KjMr ) 




if Xi > 0, and 


X = ^z!ll V ( Q n ^ 2 ) 

c_ 2d? n 5- oc L(«n -Tn/ll/d 2 ) (an+K/Mr) 

if xi < 0. Having calculated the kernel function for Im(K’) < 0, we can now let Im(A') — *■ 0. 
Then the roots become, using the branch of the square root employed previously in 
Eq. (A.37) 


e i(a„+KM r )x! (A.43) 


aj = ^ ±i^V(r„/<l) 2 - AT 2 (A.44) 

if \T n /d\ > K , and 

a± = T ^ - (r„/d) 2 ( A.45) 

if |r„/d| < AT. In Eqs. (A.44) and (A.45), a+ is to be calculated using the upper set of 
signs, and a~ using the lower set of signs. The yf~ sign means the positive square root. 
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APPENDIX B 


NUMERICAL COMPUTATION OF THE NORMAL MODES 

IN AN ANNULAR DUCT 


As discussed in Chapter 4, the normal modes in an annular duct axe functions of the 


form 


) — AJmift mn T ) 4" BYjn (^mn 7*)? (B.l) 

where J m (-) and Y m (-) are Bessel functions of the first and second kinds, and K mn ’s are 
the roots of the following transcendental equation: 


= 0 . 


J'm(Kr H ) Y^(kt h ) 

J'm{ Kr D ) Y !n( Kr D ) 

The roots K mn r D of Eq. (B.l) are found by first estimating the roots as follows: 


(B-2) 


& mn 


m, if n = 1; 

n— 1 4” if n 1- 


(B.3) 


This estimate is refined by incrementing the estimated value of Km n r D by 7r/10 until the 
determinant in Eq. (B.2) changes sign. The step size is then halved and changed in 
sign. This process continues until the absolute value of the determinant is reduced to a 
preassigned value. In the code, this is done utilizing the non-dimensional versions of the 
above equations. 


Once the eigenvalue K mn has been computed, the constants A and B are assigned one 
of the following two sets of values: 


(A = l 



{- 

J-rn ( K mn r D ) 
Ym{ K mn r b) 

or < 

< 


B = 1 


(B.4) 


Of these two sets of values, the one for which ( A 2 + B 2 ) is the smaller value is chosen. If 
(A 2 + B 2 ) is the same for both, then the second set is picked. The desired normalization, 


namely, 
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(B-5) 


is obtained by computing the value of the integral on the left-hand side of the equation 
above, using the formula 


J r)rdr = ^r 2 - ^ V'mfcrrm r) 


= c. 


(B.6) 


H 


The constants A and B are then divided by 


yj2C/(rZ - r„2) (B.7) 

to give the normalization required by Eq. (B.5). 

Note that in terms of non-dimensional variables, Eqs. (B.5) and (B.6) can be written, 
respectively, as 

l 

J i? m {Xm n x) xdx = ^(1 -<rl) (B.8) 

°T 

and 



a- 


The non-dimensional parameters in Eqs. (B.8) and (B.9) are the same as the ones originally 
used in Chapter 5. It is easily seen that the factor by which A and B are to be divided, 
can be written explicitly as 


(1 - myxj n )^(x m „) - (<7? - 


(B.io) 


- 49 - 



APPENDIX C 


SOUND POWER FLUX 

The sound power flux through any cross-section of the duct is given by (Ref. 19) 

Power = J <I ax > dA, (C.l) 

where I ^ is the axial component of the acoustic energy flux vector, < > is the acoustic 

intensity, which is defined as 


7T /Bn 


< lax > 


m f 

~ 2tt J 


fax d£, 


-tt/BO 


(C.2) 


and Ap is the cross-sectional area of the duct. In a duct containing a fluid flowing at a 
uniform axial velocity U, I ax is, in turn, given by 


Ia* = (-+Uu) ( POU + pU), (C.3) 

P o 

where po is the nominal fluid density and p, p, and u are the instantaneous acoustic 
perturbation pressure, density and axial velocity, respectively. The quantities p and u are 
real. Although not explicitly stated, all acoustic variables are functions of both space and 
time. 

The acoustic pressure and density perturbations are proportional to one another (p = 
cj}/j>), so Eq. (C.3) can be rewritten in terms of the pressure and axial velocity as 

lax = (1 + M 2 )pu+ pp + poCqMuu, (C.4) 

Poco 

where M = U /cq is the nominal axial flow Mach number. For p and u real, we can write 
pu = pu* and uu = uu* , where the superscript * denote the complex conjugate. Hence, 
Eq. (C.4) can be rewritten as 

« M 

I ax = (1 + M 2 )pu* -\ pp* + pqCqMuu*. (C.5) 

P OCO 
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For functions of £, we reverse the Fourier series sign convention used in Eqs. (3.5) and 
(3.6) for the spatial variables. Then, p can be written as the Fourier series 

OO 

P (x,t)= y, ?.(>>■“, < ce ) 

s= — OO 

where 


n/BQ 

P,(f) = J p(x,t) e“ BM dt. (C.7) 

-tt/BQ 

There is a corresponding definition for u(x,t) and u s (x ). The axial acoustic energy flux, 
lax, can now be written as a double summation 


OO OO 


M 


Lx = + M 2 )p a (x) u* {x) + —— p 3 (x)p* (x) 

— ' pqcq 

s=— oo ix=—oo u 

+ pqcqM u,(x)u1(x) 


(C.8) 




Both p a (x) and u^x) can be written as sums of modal pressure patterns as in Eq. 
(4.21). Thus 


OO OO 


V,(x) = Y E P”™ r ) 


(0.9) 


m= — oo n=l 


oo oo 


u^x) = £ £ ifc(« W r) (C.10) 

A:=— oo 

The modal coefficients p mns and u mna are not independent, but are related through the 
axial momentum equation, 


?}t + u ^L + l^L = o. 

dt dxi po dxi 


(C.ll) 


Let the pressure and axial velocity be given, for each harmonic, by the real parts, respec- 
tively, of p 3 exp(— isBQt) and u s exp(-isBQt). Then Eq. (C.ll) reduces to 


du a 

—isBClUg + U — h 


1 dp s 
dxi ' po dxi 


= 0 . 


(C.12) 
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If Eqs. (C.9) and (C.10) are substituted in Eq. (C.12), the following relationship is obtained, 


_ _ 

'U'mns — TT Pmnsi 

PqU 


where 


A Tnnx — 


Kmns 


sBCl/U + ^fmns 

The axial acoustic energy flux can now be written entirely in terms of p mn5 : 


oo oo oo oo 


/<**= XX ^fc(Kwr) e i(m fc)< ^ 


m=— oo n= 1 k =— oo £=1 
oo oo 


X E E [M 2 - (1+ M 2 ) ^ + Ann, AS<J 


po(/ X 

' 5— — OO fl=—00 


x e 


— * [(TrrmS 


Pmna Pkip.- 


Using this relation and Eq. (C.2), then 


OO OO OO OO 


</»>= E E E E wiw'iAMf"- 11 * 


m= — OO TT— 1 fc= — oo £=1 


SO 


OO OO 

XX / “ (•*■ + -^ 2 ) + x4 mns 

=-oo M =-oo _^ Bn 

Pmn3 p * dt = 


t r/£fi 


2npoU z 

5=-oo M=-oo _ w/Bn 


x e 

oo oo oo oo 


XXX XI *pTn{Kmnr) ll>k(Kkir) e i(m *>* 


m= — oo n=l /c=— oo £=1 


X -U [A* 2 - (1 + M 2 ) AJ (j + Ann, AkJ 
X e - i(7m " s - 7fc<a)l1 p mn , p* kta . 

Note that in obtaining Eq. (C.16), those terms in the integrand for which /x f 
periodic in t over the interval 2ir/BQ, and thus they average out to zero and 
contribute to the integral. 


(C.13) 


(C.14) 


(C.15) 


(C.16) 


s, are 
do not 


- 52 - 



Substituting for < 1^ > in Eq. (C.l), the sound power is found by performing the 
integration over the cross-section of the duct. Utilizing the orthogonality conditions that 
exist, to eliminate product terms involving different modes, we find that 


2tt r D 

Power = 

0 r H 



< lax > r dr d<t> = 


*(r p 2 -r„ 2 ) 

PoU 


(C.17) 


oo oo oo 

X y] y] [M 2 — (1 + M 2 ) + (Annol ] |pmns | j 

m=— oo n=l 3— — oo 

where the integral over the cross-sectional area of the duct is written explicitly in terms 
of r and <f> coordinates. 

The term in the square brackets which, henceforth, we denote as can be sim- 

plified, resulting in 


^M 2 p*(sBQ/U)k 

TTXn3 fr 1 1 

mns — ~Z 7~2~ ’ 

[sBfl/ Co ± Mk m na] 

where the upper set of signs applies upstream of the stator and the lower set applies 
downstream. In terms of G mns , then, we have 


Power = 



OO OO OO 

y yy yy G mn5 \pmns\ • 

m= — oo n — 1 3= — oo 


(C.19) 


Finally, recalling that the index m is related to the summation indices s and q through 
the relation m = sB - qV, Eq. (C.19) becomes 


- j) 

PoU 


oo 

E 

oo 


E E 

n=X oc 


|Pmnj | 


Power = 


(C.20) 



APPENDIX D 


FILON’S INTEGRATION RULE 

Filon’s integration rule applies to integrals of the form 

b 

J f(x)e i9{x) dx, (D.l) 

a 

where the phase function g(x) is a “large” linear function of x, so that use of the trapezoidal 
rule or Simpson’s rule would require that the interval of integration be divided into many 
subintervals to obtain an accurate answer. Filon’s rule is obtained by assuming that the 
function f(x) can be approximated by a quadratic function, but the exponential function 
exp[ig (x)] is integrated exactly. To extend this procedure to integrals wherein g(x) is not 
a lin ear function of x, we need only approximate it as a linear function within each x 
subinterval. The most straightforward procedure is to approximate both f(x) and g(x) as 
linear functions of x. For example, to compute the following integral 

b 

1 = J /(*) e i9{x) dx, 

a 

divide the interval (a, b ) into N equal subintervals of length h 
over the nth subinterval is then 

x n +h 

In= J f(x)e i9 ^dx. (D.3) 

To compute this integral, approximate both f(x) and g(x) as linear functions: 

m = /„ + (/„+, - . (°- 4 ) 

S (*)= 9 „ + (9„ +1 -S„)(^), (D.5) 

where f n — f(x n ), and so on. Carrying out the integrations, we obtain 
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(D-2) 

= (b — a)/N. The integral 



In = Ma'fn + a/ n +i) e'», 


(D.6) 


where 


50 = 2 (5n + 5n+l)j (D-7) 

A = ^(5n -5n+l), ( D - 8 ) 


sin A / sin A 
~ 2A + %2A 2 



(D.9) 


and a* is the complex conjugate of a. When A — * 0, the integral over the nth subinterval 
becomes 


In = 2 (fn + fn+ l)j (D.10) 

which is the trapezoidal rule. 

If an even number of subintervals are used, it is also possible to approximate f(x) as 
a quadratic function (over any two neighboring subintervals) while leaving g(x) as a linear 
function of x. This integration scheme reduces to Simpson’s rule when g(x) is constant. 
Both the trapezoidal and the quadratic versions of Filon’s rule have been tried. Because the 
quadratic version did not seem to improve the convergence significantly over the simpler 
trapezoidal rule, the latter scheme is used in the V072 program. 
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APPENDIX E 


NUMERICAL SOLUTION OF THE INTEGRAL EQUATION 
FOR THE BLADE LOADING 

As discussed in Chapter 5, the elemental blade loading function f a {y/b) is the solution 
of Eq. (5.8), which with subscript s suppressed in this appendix, becomes 

e Al/i = [' Kci^) Hr) <*(£)• (E.1) 

Jy/b =- 1 b b b 

This integral equation contains two basic difficulties. They are: 

(a) the solution f(z/b) has a singularity of the type ( z + 6)“5 
at the leading edge (z = —b), and 

(b) the kernel function K c {y/b ) contains both a Cauchy 
singularity (1 /y) and a logarithmic singularity at y = 0. 

Difficulty (a) is circumvented by introducing the chordwise transformations z/b = cos 6 
and y/b = cos ip. Equation (E.l) then becomes 

K c ( cosO — cos ip) F(ip) dip , (E-2) 

where F(ip) = /(cos ip) sin ip. Whereas f (cosip) is singular at ip = 7r, F(ip) is not. Thus 
the integral equation is solved for F rather than for / itself. Difficulty (b) is overcome by 
special treatment of the logarithmically singular part of the kernel. 

To solve Eq. (E.2), we apply the method of collocation. That is, the integral on the 
right-hand side of Eq. (E.2) is required to equal the forcing function on the left-hand side 
at the N points 


Ak cos 0 


=/ 

Jo 


Om = (m - 1/2)7t; m = 1, • • • , TV. 
The equations to be solved are then 

e ik cos e m _ j(Q m )- m= 1 ,..., N, 


(E.3) 


(E.4) 
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where 


1(0) = J K c (cos9 — cosxp) F(xp) dip. (E.5) 

o 

In computing 1(9) by evaluating the right-hand side of Eq. (E.5), the points 

iPe=£ir/N, £ = 0,1, - • • ,N, (E-6) 

are used.* Having the collocation and integration points interlaced has the dual effect of 
avoiding the point xp — 9, where K c is singular, and of ensuring that the solution obtained 
satisfies the Kutta condition, which requires that F = 0 at the trailing edge. 

The logarithmic portion of the kernel function must be isolated for special treatment. 
Thus, 


and 


where 


KJy/b) = Ki(y/b) + G(y/6) log |y/6| 


m=hw+hm, 


(E.7) 


(E.8) 


7T 

l(*) = / K\ (cos 6 — cos rp) F(ip) dip, 


(E.9) 


7T 

Il(9) = J G(cos0 — cosTp)\og\cos0 — cosxp\ F(xp) dip. (E.10) 


The trapezoidal rule suffices to compute h(9), 


N 

h(9) = Y, jj B e K x (cos9 - cos fo) F(^), 


e=o 


(E.11) 


* Because F(0) = 0, use of N + 1 integration points as called for in Eq. (E.6) introduces 
only N unknowns F(£ir/N), for £ = 1, • • • , N. 
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where Bo = Bn = 1/2, while Be = 1 otherwise; but a special integration rule, devised by 
Whitehead (Refs. 20, 21) is required for the logarithmic portion of the kernel. Whitehead’s 
integration rule is given as 

N 

I L (0) = ^2 ^B/G(cos0-cos^)log|cos0-cosV>|S/(0), (E.12) 

e=o 1 

where the B/s are the same weighting functions as defined above, and 

Sg(0) = -|log2 + 2^ ^B s cos(s^) cos(s£7r/iV)|. (E.13) 

Note that because of the comment in the footnote on page 57, the summations in Eqs. 
(E.ll) and (E.12) actually begin at t = 1 rather than l = 0. Hence by combining I\(0) 
and Il(0), the integral equation is reduced to the following set of algebraic equations for 
the unknowns F(£n/N), 

N 

e ik cos{(m— l/2)7r/iV = ^ A^F^/N), TO = 1, • • • , N, (E.14) 

1 = 1 

where 

Am = i B e I.K c (x) + G(i) [S,(« m ) - log |x|] ) (E.15) 

and 

x = cos(0 m ) — cos (ipe)- (E.16) 

Any standard matrix package which handles equations with complex coefficients can be 
used to solve these equations. 
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